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We investigate how a Type IV future singularity can be included in the cosmological evolution 
of a well-known exponential model of inflation. In order to achieve this we use a two scalar field 
model, in the context of which the incorporation of the Type IV singularity can be consistently 
done. In the context of the exponential model we study, when a Type IV singularity is included 
in the evolution, an instability occurs in the slow-roll parameters, and in particular on the second 
slow-roll parameter. Particularly, if we abandon the slow-roll condition for both the scalars we shall 
use, then the most consistent description of the dynamics of the inflationary era is provided by the 
Hubble slow-roll parameters ch and i]h- Then, the second Hubble slow-roll parameter rjH, which 
measures the duration of the inflationary era, becomes singular at the point where the Type IV 
singularity is chosen to occur, while the Hubble slow-roll parameter ch is regular there. Therefore, 
this infinite singularity indicates that the occurrence of the finite-time singularity is responsible for 
the instability in the scalar field model we study. This sort of instability has it’s imprint on the 
dynamical system that can be constructed from the cosmological equations, with the dynamical 
system being unstable. Also the late-time evolution of the two scalar field system is studied, and in 
the context of the theoretical framework we use, late-time and early-time acceleration are described 
in a unified way. In addition, the instability due to the singularity mechanism we propose, is 
discussed in the context of other inflationary scalar potentials. Finally, we discuss the implications 
of such a singularity in the Hubble slow-roll parameters and we also critically discuss qualitatively, 
what implications could this effect have on the graceful exit problem of the exponential model. 

PACS numbers: 04.50.Kd, 95.36.-l-x, 98.80.-k, 98.80.Cq,11.25.-w 


I. INTRODUCTION 

sin- 
link 
two 

alternative classification schemes, the one that uses the scale factor as a rule for the classification and another 
classification scheme that uses the scale factor, the effective energy density, the effective pressure and in some cases 
the Hubble rate in order to determine the type of the singularity Si- The most interesting singularities which 
are typical for quintessence phantom dark energy evolution are finite time future cosmological singularities and they 
were firstly developed in [I| and in @,0 respectively. The most severe, and rather unwanted types of singularities 
in every theoretical physics framework, are the crushing type singularities, like the initial singularity or for example 
the Big Rip d, 0. For these singularities, the Hawking-Penrose strong energy theorems [1| are violated and therefore 
geodesics incompleteness occurs. These singularities are possibly indicators of the inability of the theoretical frame¬ 
work that predicts them, to fully describe the physical evolution consistently at the spacetime points that these occur. 
So probably these can be viewed as possible indicators of the compelling need of a more fundamental theory for the 
consistent description of the physical phenomena near these spacetime points. For some informative studies on this 
kind of singularities, see |9l-Hll| and references therein. Also for insightful studies on the avoidance of the initial singu¬ 
larity see There exist however other types of singularities that can be characterized as ’’mild” singularities, 

with mild indicating the fact that no geodesics incompleteness occurs at these spacetime singularities and also the 
Hawking-Penrose theorems, and even the weak energy theorems are not violated. One characteristic example of this 
type of singularities is the Type IV singularity, which can be viewed as a link between the classical and the quantnm 
description (whatever it may be) of the theory under study. Some properties and also phenomenological-observational 


One of the unavoidable features of many cosmological theoretical models is the appearance of cosmological 
gularities, which either signal the need for an alternative physical description or maybe can be viewed as the 
between a classical and a quantum cosmological theory. In cosmology singularities are considered according to 


2 


implications of Type IV singularities as applied mainly to inflation, were recently studied in IldMlTII . 

Particularly, in Ref. 0, using standard reconstruction schemes for scalar-tensor theories |l8l - l20l| . we incorporated 
a Type IV singularity in a model of inflation, by using two scalar fields. As we demonstrated, the second scalar 
field has no effect in the early-time dynamical evolution of the physical system and contributes only at late-times. 
The purpose of this paper is two fold: Firstly we shall demonstrate how we can consistently incorporate a Type IV 
singularity in a quite well known exponential scalar model of inflation [HHil of the form V{ip) = Vbe-^"^, by using 
two scalar fields. Secondly, we shall demonstrate that a singularity occurs in the Hubble slow-roll parameters of the 
exponential scalar model with potential V{}p) = Note that according to the Planck data the model we 

shall study is consistent to some extent to the measured observational data, but there is no exit from inflation for this 
model, so the instability we found can maybe have an impact on the graceful exit problem of the exponential model. 

Using two scalar fields, we shall demonstrate that the theory has a source of instability for this scalar model. 
Particularly, the main reason behind this instability is the Type IV singularity. As we shall evince, owing to the Type 
IV singularity, the theory becomes unstable because the slow-roll parameters and specifically the second slow-roll 
parameter develops an infinite singularity at the Type IV singularity point. Then, by suitably choosing the time at 
which the Type IV singularity occurs, we can render the system unstable at any point we might wish. Note that 
the second slow-roll index measures how long the inflationary era is, while the first slow-roll parameter indicates if 
inflation occurs. In the context of the two scalar field model we will work on, the second scalar has a negligible effect 
at the early time evolution, apart from the instability that it causes via the Type IV singularity, and the early time 
evolution is solely governed by the canonical scalar field with scalar potential the exponential potential of the form 
V{tf) = In addition, the second scalar will eventually control the late-time era, although this can be strongly 

model dependent. One of the important assumptions we shall make is that the canonical scalar field corresponding 
to the exponential potential V{ip) = , does not satisfy the slow-roll conditions, in contrast to the case studied 

in our previous study [l^ . in which the slow-roll conditions were determined by the potential. Therefore, in that 
case, the singularity did not modify or affect the slow-roll conditions, unless the time that the singularity occurs was 
chosen before the inflation ending time. The definition of the slow-roll parameter using the potential was described 
in detail in l?7| an d was firstly given in [2^. In contrast, the slow-roll condition can involve the Hubble rate, as was 
firstly done in [29|, and we adopt this approach in this paper. As we shall demonstrate, the singularity in the second 
slow-roll index is due to the appearance of the second derivative of the Hubble rate in it’s functional form. With 
regards to the dynamical system instability caused by a singularity there, we have to note that this sort of instability 
is somehow different from the infinite instabilities occurring in the effective equation of state of single scalar field 
models when the phantom divide is crossed. These instabilities were firstly pointed out in [30l| and were studied in 
detail in [mill. Finally, we demonstrate that in the context of the two scalars model we shall use, it is possible to 
unify the early-time with late-time acceleration. 

This paper is organized as follows: In section H we present in brief the classification of finite time singularities and 
also the geometric background conventions we shall use throughout the rest of the paper. Moreover, we describe the 
basic features of the exponential scalar model, including the observational predictions that it implies. In addition, we 
present the two scalar field model and study it’s dynamical evolution by appropriately constraining the second scalar 
field. In section HI we present in detail the two sources of instability and also we qualitatively discuss the potential 
implications that the instability could have on the graceful exit problem of the exponential potential. Finally, the 
late-time dynamics of the two scalars model is studied in section IV and in section V we present the implications of 
the Type IV singularities to other scalar models. The concluding remarks along with a discussion of the results follow 
at the end of the paper. 


II. EXPONENTIAL TWO SCALAR MODEL OF INFLATION AND TYPE IV SINGULAR EVOLUTION 

As we already mentioned, the detailed classification of finite time singularities was firstly developed in Refs, 
and we now present in brief the essential features of this classification scheme. 

• Type I (“The Big Rip Singularity”): This singularity occurs when the cosmic time approaches t —>■ ts, where 
the following physical quantities behave as: the scale factor a, the effective energy density peff and finally the 
effective pressure Peff diverge, that is, a —>■ oo, pes oo, and |peff| —>■ oo. This singularity is of crushing type, 
and more details can be found in Refs. 

• Type H (“The Sudden Singularity”): This singularity occurs when the cosmic time approaches t —>■ fy, where 
only the scale factor a and the effective energy density peS tend to a finite value, that is, o —>■ Os, PeS —>■ Ps, but 
the effective pressure diverges, |peff| —>■ oo. For details on this singularity, consult Refs. Si- 

• Type HI : This singularity occurs when the cosmic time approaches f —>■ tg, where only the scale factor tends 
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to a finite value, that is, a —> a^, but the rest two physical quantities, that is, the effective pressure and the 
effective energy density, both diverge, that is, |peff| oo and Peff c»- 

• Type IV : This singularity is not of the crushing type and it is the most ’’mild” among the four types of finite 
time cosmological singularities. In this case, as t > tg, all the aforementioned quantities tend to finite values, 
that is, a —>■ Os, PeS —>■ Ps, |peff| —>■ Ps, and also the Hubble rate and it’s first derivative are finite. But the 
second or higher derivatives of the Hubble rate diverge as t —>■ ts ■ For details on this type of singularity see Ref. 

9 . 

Having presented the essentials of the finite time cosmological singularities classification, it is worth presenting 
also in brief the geometrical conventions we shall use in the rest of the paper. We shall assume that the background 
spacetime metric is a flat Friedmann-Robertson-Walker (FRW), with line element of the form, 

ds^ =—dt^ + a{t)^ ^ {dx''Y , ( 1 ) 

1=1.2,3 


where the parameter a{t) stands for the scale factor. Also, the effective energy density peff and the effective pressure 
of the matter fluids Peff for the FRW of the form 0 are related to the Hubble rate in the following way, 

PeS = ^H\ peff = -4(2^ + 3Ff") ■ (2) 

Originally, the finite-time future singularities were studied after dark energy era. However, recently it was realized 
that they are relevant and some of them maybe quite realistic also during or just after inflation. In our previous 
studies |14l - ll7j| we incorporated a Type IV singularity in the cosmological (mainly inflationary) evolution of some 
scalar-tensor models using one scalar [T^ or two scalars [Il[i3. But in the case of the nearly inflation potential 
which we studied in [l^ . the use of two scalar fields was compelling, for reasons we explained in detail in [l^ . The 
same reasoning applies in the case we study in this paper too, therefore this incorporation can be realized by using 
two scalar fields (details on the reasoning why to use two scalar fields are presented in the appendix). Before going 
into the details of our model, we shall present in brief the exponential model of inflation [2ll - [^ and we also describe 
when this model can be partially compatible with Planck data. Also we shall demonstrate why it is impossible to 
have a graceful exit from the inflationary era, regardless if we assume slow-roll or not. 

The exponential model we shall be interested in, is a power-law model of inflation in which, the action of the 
canonical scalar field is the following, 

S = I - ^(<P)I : (3) 

with the scalar potential being of the following form, 

V{p) = , (4) 

and with = 8ttG, where G is Newton’s constant. The exponential potential of Eq. o, has been used in many 
cosmological contexts and for an incomplete list on this see [2ll - l^ . As we shall present in detail later on, there exist 
two alternative descriptions for the slow-roll parameters, one that takes into account the potential [13, and the 
other one takes into account the Hubble rate only |23iJ3. By taking into account only the scalar potential, the 
slow-roll parameters are defined as follows, [13, [H, l32l l33j|. 


1 (v'ip)\ 

2^2 Vip) ) 


k2 \V{p) j 


(5) 


The observational indices of inflation can be defined in terms of these slow-roll parameters, and we shall be interested 
in the spectral index of primordial curvature perturbations Us and the tensor-to-scalar ration r. These are given in 
terms of the slow-roll parameters in the following way [13 . [33l| , 

r = 16e , Os ~ 16e?7 — 24e^ — 2^^ . (6) 



ris ~ 1 — 6e -I- 2ry 
For the potential (|3]), the slow-roll indices read. 


1 

e = - 

q 


( 7 ) 
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which are called in the literature ’’potential slow-roll parameters” [27|, so we adopt this terminology too. The 
corresponding observational indices are easily calculated, and these read, 


ris ~ 1 - 


16 


( 8 ) 


However, in this paper we shall not assume that the slow-roll approximation holds true, and therefore we need an 
alternative description for the corresponding slow-roll parameters that measure the inflationary dynamics. As was 
discussed in , the slow-roll parameters of Eq. 0, when these are small, provide a necessary consistency condition 
for inflation to occur, but not a sufficient one. A much more superior definition of the slow-roh parameters 27[ that 
encompass all the inflationary dynamics, are the Hubble slow-roll parameters eu and r]H [271 l29l|. which are defined 
in terms of the Hubble parameter as follows, US, [11, 




2 


, VH = 


2_EM 


(9) 


where again we used the terminology of [27| . It is more convenient for the purposes of our analysis to express the 
Hubble slow-roll parameters in terms of time and following , these are given by. 


Cff = 3 


H(^) + f 


■qu = -3 




ZHif 


Using the following equations of motion for the canonical scalar field i^, 

Sttk^ / 


+ V{‘p) , (p + iHiP + V'{ip)=Q, 


^ = 3 V2 

the Hubble slow-roll parameters (US can be written in the following way. 


H 


H 


en = 


'i?2’ 2HH 


( 10 ) 


( 11 ) 


( 12 ) 


In the rest of the paper we will make use of the Hubble slow-roll parameters defined in Eq. m, which are suitable 
for our analysis. Note that in our previous work, related to the inflation potential, we made use of the potential 
slow-roll parameters, but in that case, the graceful exit was guaranteed by the breaking of the slow-roll conditions, a 
very well known mechanism [s^. For the potential (HJ, the Hubble slow-roll indices read. 


1 1 

e = V ^ 

9 9 

and the spectral index of primordial curvature perturbations reads 

n., = 1 - 


9-1 


(13) 


(14) 


The latest Planck data indicate that the values of the aforementioned spectral indices are constrained as follows, 


n* = 0.9644 ±0.0049, r<0.10, (15) 

so the scalar-to-tensor ratio is already excluded (see also HI), but the index Ug can be compatible to the 2015 Planck 
data jl^, if the parameter q takes the following values, 

0<9<101, (16) 

so we choose q = 100 for our analysis. Before we proceed, let us discuss the issue of graceful exit from inflation for 
the model (jll). By looking at the potential slow-roll parameters and also at the Hubble slow-roll parameters, in both 
cases it is not possible to have graceful exit since these never become of order one, and also in the context of a single 
scalar field model, there is no alternative mechanism to trigger the graceful exit from inflation. We have to mention 
that in the literature there are two very popular mechanisms that can trigger the exit from inflation, with the first 
being related to tachyonic instabilities |24| and in the context of the second mechanism, the graceful exit is triggered 
by a trace anomaly |35l |. 
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In this paper, we shall present a mechanism that renders the system dynamically unstable, with the instability being 
triggered by the Type IV singularity. This mechanism requires the use of two scalar fields, as we now demonstrate. 
In the appendix we provide the details on why it is very difficult to incorporate the Type IV singularity to the single 
scalar field exponential model with potential (|H). 

In order to incorporate the Type IV singularity in the cosmological model described partially by a canonical scalar 
field with scalar potential that of Eq. o, we shall use a two scalar field scalar-tensor model. In principle, in the 
context of two or more scalar field models, there are much more options for a successful description-realization of 
various cosmologies. In addition, the two scalar field models can remedy certain inconsistencies that arise in the single 
scalar field models, mainly having to do with the effective equation of state (EoS) of the scalar field when the phantom 
divide crossing occurs [1^ [l^ . Among the vast possibilities we can choose for the double scalar field model, we 

shall make a convenient choice that proves to be valuable from a phenomenological point of view. 

Consider the following two scalar field action, 

S = J ~ \v{x)d^,xd^X ■ ( 17 ) 

where as usual, uj{4>) stands for the kinetic function of the non-canonical scalar field (j), while ??(x) denotes the kinetic 
function of the second non-canonical scalar field y. Of course it is conceivable that when one of the functions 
or ri{x) is negative, then the associated to this function scalar field becomes a phantom field, but as we shall see, 
by appropriately choosing the initial conditions of the scalar field y, the possibility of a phantom scalar field and 
consequently of partially phantom inflation is avoided. Notice however that this issue is strongly parameter and 
model dependent. Assuming that the scalar fields depend only on the cosmic time t, then for the FRW cosmological 
background o, the FRW equations for the action (13 read, 

+ v{x)id =-\h , + h) . (18) 

AC AC \ / 

By choosing the generalized scalar potential V{(j),x) and the functions uj{(j)), r]{x) to satisfy the following relations, 

oj{t)+r]{t) =+ f{t)) , (19) 

AC AC 

we can find an explicit and quite general solution of Eqs. dni), which is of the form, 

0 = y = t, H = f{t). (20) 

Practically, what we i^t described is the scalar reconstruction method, and for further details and applications, the 
reader is referred to |18l| . Since this method provides us with much freedom for choosing the field kinetic functions, 
we choose these as follows. 


a^(<(>) = {/'(<(>) - \/a((^)2 -p/'(0)2| , r?(y) = -^V<^{x)^ + fix)^ ■ (21) 

where the function a(x) is arbitrarily chosen and it’s explicit form will be specified shortly. The scalar potential can 
be specified in terms of an auxiliary function /(<)), y), which is defined to be equal to, 

Ww) = -y (^J d 4 )u:{ 4 >) + J dy??(y)^ . (22) 

and is chosen to satisfy the following, 

(23) 

Then, the scalar potential can be expressed in terms of the auxiliary function /(</>, y), as follows. 


x) 


1 


^3/((/),y)^ 


X) 

dcj) 


x) \ 

9x j 


(24) 


Notice that, the constraint ( 1 ^ 51 ) fixes the values of the integration constants that arise in Eq. ( 1 ^^ . The two non- 
canonical scalar fields (j) and y satisfy the following equations of motion, 

0 = + l-uj '® + lv'ix)x^ + ^Hr]{x)x + 

2 d(f> 2 dy 


( 25 ) 
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In order to consistently incorporate a finite time cosmological singularity in the two scalar field cosmological evolution, 
the Hubble parameter is assumed to have the following general form, 

Hit) = h{t) + f2it)it,-tr, ( 26 ) 

where the parameter a actually determines the type of the finite time cosmological singularity, and it will be assumed 
to take non-integer values of the form. 


a = 


n 

2m + 1 ’ 


( 27 ) 


with m and n integer numbers, the values of which will be determined shortly. According to the classification of finite 
time cosmological singularities we presented previously, the values of a classify the cosmological singularities of the 
evolution model (j26ll in the following way, 

• q; < — 1 corresponds to the Type I singularity. 

• — 1 < q; < 0 corresponds to Type III singularity. 

• 0 < q; < 1 corresponds to Type II singularity. 

• q; > 1 corresponds to Type IV singularity. 

Since we are interested in the Type IV singularity, for the rest of this paper it is assumed that a > 1, so the integers 
m and n appearing in (1271) . can be chosen in such a way so that a > 1. For the purposes of this paper, a is assumed 
to be equal to a = 4/3, for reasons that will become obvious soon. 

A convenient choice for the function a{x) which appears in Eq. (EH), is the following, 


a{x) =+ af 2 {x){ts-x)°‘ ( 28 ) 

and it is conceivable that the variable x stands for the scalar fields (j) or x- By using this form of a(x) and substituting 
it to Eq. EH, the kinetic functions of the non-canonical scalar fields become significantly simplified and are equal to, 

w((/) = , "nix) = -^ (fLix) {ts - x)“ + af2{x) (ts - X)“"^) , ( 29 ) 

AC K, \ / 


In addition, by using ( 1 ^ . the auxiliary function f{(j),x) of Eq. ( 1 ^ . becomes. 


x) = /!(</) + f 2 (x) (ts - x)“ , (30) 

and consequently, the two scalar field potential V (</>, x) becomes, 

^ ifiW + f 2 {x) {ts - xff + + fLix) {ts - xT + af 2 {x) {ts - xT~^) ■ ( 31 ) 

Equipped with Eqs. EH and EH, we can easily realize a singular evolution related to the scalar potential (|1])- 
Specifically, we choose the Hubble rate (1261) to be equal to, 

H{t) = + C4 (t« - tf , (32) 

C2 + Cst 

with the parameters c^, i = 1, ...,4, being for the moment arbitrary constants to be specified soon. Using Eq. (I32L 
the function a{x) appearing in Eq. ( 1 ^ 51 ) becomes equal to. 


y{x) = \ cl{ts - x) 2-l-2aQ,2 _ | _ 


ClC3 


(C2 -I- C3x)2 

and also the kinetic functions a;((/)) and r?(x) of Eq. (1^ become equal to 


- C4(<s - x) i+“c 


;(</) = 


2ciC3 


K^{c2 + C3(/)2 


vix) = - 


2c20L{ts - x) 


— l+a 


( 33 ) 


( 34 ) 
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Notice that in principle, the scalar field x can become a phantom scalar, and indeed this can be the case ^ when 
a = 4/3. However, as we will demonstrate, the scalar field x turns to phantom much more earlier than the infiationary 
era. In addition, in view of the choice (I33L the auxiliary function /(/>, x) takes the following form, 


/(</>, x) = 


Cl 


C2 + C3<t) 


+ CA{ts - XY 


and therefore, the scalar potential of Eq. m reads, 


mx) = 


3cf 


C 1 C 3 


C4a(ts - x) , 6ciC4(ts - x)“ , 3cl(ts - x) 


2a 


k2(c2+C3(/))2 k2(c2+C3(/)2 


k 2 (c 2 + C 3 (/)) 


(35) 


(36) 


In order to obtain the scalar potential of Eq. o, we will transform the scalar field </> to a canonical scalar field (p, by 
making the transformation 


rv 

ip= d(j)\J ijj{(j)) 

with being of the form given in Eq. (IMl) . which yields. 


C2 + C3(/) = e V 2<=i 

By using the canonical scalar field the action of Eq. (EZD takes the following form. 


- r 1 


1 




S' = y 

with the scalar potential V{p,x) being equal to 


1 /2c2a(x - ts) 


— 1+a 


df,xd^X-y{v’,x) (■ ■ 


(37) 


(38) 


(39) 


= fif - fiSj) - V 

\ K,^ J 


— l+a 


6 ciC 4 (ts - x)“e 


+ 


- xY 


(40) 


where, since a = 4/3, we used the following relation in Eq. (15^ . 

(i.-x)-'+“ = -(x-ts)-'+“ (41) 

It is obvious that the potential of Eq. (HOI) contains the scalar potential of Eq. plus the potential of the second 
scalar field x- In order for the exponential potential of Eq. (l40l) to be identical to the scalar potential of Eq. o , we 
must require that. 


Eo = 


3cf 


C 1 C 3 
2 ’ 


Cl 

— = 9 - 

C3 


(42) 


At this point we shall make a critical assumption for the dynamical evolution of the system of the two scalar fields p 
and X, described by the potential (l40l) and for a Hubble rate given in (1321) . Particularly, we require the following. 


• The canonical scalar field p does not satisfy the slow-roll requirements. 

• The non-canonical scalar field x does not satisfy the slow-roll requirements too. 

• The non-canonical scalar x is assumed to have very small values at early times. 

• The parameters c^, i = 1,.., 4 are chosen in such a way, so that at early times the contribution of the scalar field 
X to the potential V{p, x) given in Eq. (HOI) is negligible, and also the kinetic fnnction of the scalar field x is 
also negligible at early times. 


^ It is easy to see this, since for a = 4/3, the function ri{x) is equal to, ??(x) = -’ which can be phantom 

when X < is- For the values of the parameters that we will choose however, this never occurs 
















• The Type IV singularity, which occurs at the time tg, is chosen to appear at an early time 10 ^®sec, where it is 
usually assumed that inflation ends [s^. 

By taking into account the requirements of the list above and combine these with the constraints (l42l) . and also that 
we assumed q = 100, we choose the parameters Ci, i = 1,4 as follows, 

Cl = 10^^ C 2 = 10-2°, c3 = 10^°, C 4 = 10-‘‘° (43) 

Also we shall take tg, to be equal to tg = 10-°°sec, which is the time that inflation ends in most theories. Of course 
inflation does not end in the present model for the moment, but as we describe soon, there is a mechanism that can 
potentially create a dynamical instability and this instability could be an indicator of graceful exit. 

For the choice of the parameters (l43l) . the Hubble rate is positive, as it can be seen in Fig. O- 



FIG. 1: The Hubble rate H{t) = + C 4 {tg — t)°‘ as a function of the cosmic time with ci = 10 ^^^ cj = 10 , C 3 = 10^°, 

C 4 = 10—^°, ts = 10“^®sec and a = 4/3. 


Before continuing our analysis we need to discuss the issue of the physical units we shall adopt. For the purposes of 
our analysis, we shall adopt unit system for which h = c = 1, where c is the speed of light. We shall choose to express 
every physical quantity in seconds. So the length I is measured in seconds (c = 1), time in seconds, mass is measured 
in sec”^ (since I = h/mc = 1/m, when c = 1). Therefore, the canonical scalar field is measured in sec”^, the Hubble 
rate is measured in sec”^, while the potential in sec”'^. Finally, the non-canonical scalar field is dimensionless, since 
^2 is measured in sec2 and the derivative in sec”^. In the following, we shall use these dimensions in our analysis and 
it is conceivable that all parameters can be expressed in these units. 

For the values of the parameters (1431) . and by choosing a significantly small initial condition for the scalar field x 
at early times (earlier than tg = 10-°°sec) the contribution of the scalar field to the potential V{(p,x) of Eq. (I40L 
is negligible. Also the kinetic term of the scalar field x is also negligible at early times. We assume that the initial 
conditions for the scalar field x at ^ = 0 is, 

x(0) = 10-25, x'(0) = 10, (44) 

for reasons which will become clear when we discuss about the instability issue. Having all these considerations into 
mind, we conclude that at early times, the scalar field action is governed by the canonical scalar field (^, and is 
approximately equal to, 

S = J d'^xy/^ { - F((p) I , (45) 


with V{ip) given in Eq. (jT]), since the term _|_2£i£±LiiAi_£-and also the rest x-dependent terms in Eq. (HOI) , are 

negligible. In addition, as the cosmic time t grows larger, the contribution of the second scalar field x will eventually 
dominate the dynamics of the model, so that at late times, the scalar action is approximately equal to. 




1 / 2 c 2 a{tg - x) 

2 V ^2 


duXd'"x-Vxi‘P,x) 


(46) 


with Vy^{ip,x) being this time equal to. 


Pr t \ Cia{tg - x) , 6ciC4(ts - x)“e 
^xW,x) ---0-+--2- 


+ 


3cl(t. - x)"“ 


( 47 ) 
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Depending on the dynamics and on which time we choose inflation to end, it is possible that the (^-dependent 
exponential in (1471) can be either equal to one, or it can be approximately equal to zero. The qualitative picture in 
both cases is the same, but the quantitative picture is different, with the scalar field y becoming dominant much more 
earlier in the case that the (/^-dependent potential term is included in the analysis. So it is important to study both 
cases. 

Before proceeding to the two mechanisms that can generate the instability of the system under study, we study 
numerically the dynamical evolution of the scalar field y, for the choice of parameters (|43l) . in order to verify that 
our claims can be realized in the model of the two scalars. We will show that the scalar field y provides a negligible 
contribution to the two scalar fields model action at early times, but can dominate the late-time evolution. 

We intend to study in detail the dynamical evolution of the scalar field y, by solving numerically the following 
equation of motion, 


vix)x + \v'{x)x^ + ^Hr]{x)x + = 0 (48) 

2 dx 

which governs the non slow-roll evolution of the scalar field y. Recall that the potential V{(p, y) is approximately 
equal to. 


Viv’,x) - 


C4a{ts - y) GciC4{ts - y)“e 3cl{ts - y)^“ 

o o o 


(49) 


6c c (t — )°‘e~ y 

if the scalar field tp is assumed to have quite large values at early times, so the term -|_ 2 £i£ib. 2 _^ij£-is neglected. 

We assume that the scalar field y satisfies the initial condition (HTl) . in conjunction with y'(0) ~ fO, and that the 
parameters are chosen as in Eq. (H51) . Also note that the present time in seconds is approximately tp ~ 4.25 x fO^^sec 




FIG. 2: The time dependence of the scalar field x(t), with ci = 10^^, C2 = 10 , cs = 10^°, C4 = 10 ts = 10 ®®sec, a = 4/3 

and for the initial conditions y(0) ~ 10“^®, X^(0) ~ 10. 

and also that = SttG = 2.0944 x 10“^®GeV“^, with G denoting Newton’s constant. Using these, in Figs. [2]and|3l 
we present the time dependence of the scalar field y. As it is obvious, the scalar field y is very small at early time, 
and remains small for a wide range of time, until present time where it grows large. After the present time era, it 
evolves larger and dominates the evolution. This qualitatively appealing behavior was also pointed out in [T^ . and 
we briefly address this issue in a later section. It is worth checking whether the inclusion of the (^-dependent potential 



t (sec) 


FIG. 3: The scalar field’s x(t) evolution as a function of cosmic time t with ci = 10^^, C 2 = 10 , cs = 10^°, C 4 = 10 

ta = 10“®® sec, a = 4/3 and for the initial conditions x(0) ~ 10“®®, X^(0) ~ 10. 
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term 6eic4(ts x) 


in the scalar field x potential resnlts to a different behavior. In the most extreme case, the 


(/3-dependent term will be of the form ^ -which occurs if ip is approximately equal to zero after inflation ends 

(assuming that our proposed mechanism ends inflation eventually). However, even including this term, the qualitative 
behavior remains the same, as it can be seen in Fig. |4j but the scalar field x dominates to a great extent the evolution 
much more earlier in comparison to the other case. The same analysis can be performed for the kinetic function r]ix)i 



FIG. 4: The time dependence of the scalar field x(t), with ci = C2 = 10 cs = 10^°, C4 = 10 , ts = 10 ®®sec,a = 4/3 

and for the initial conditions x(0) ~ 10“^®, X^(0) — 10, including the (p-dependent term. 


which as we claimed previously, can be neglected at early times. Indeed, in Table HI we have calculated the values of 
the kinetic function ri{x{i))i for various time moments and also for the cases that the (/^-dependent potential term is 
included (indicated by ?7u,(x(t))) or not (indicated by ry„(x(t))). As it can be seen from Table (|T|, the behavior of the 


Time 

1 

O 

t ~ 10"^® 

t ~ 10“® 

t ~ 10“'‘ 

ViXnit)) 

3.48 X 10"^® 

2.35 X 10”^° 

H- ' 
H- ' 

o 

1 

-q 

8.8 X 10"'^ 


3.47 X 10"^® 

2.2 X 10”^° 

8.8 X 10“ 

9.2 X 10^® 


TABLE I: for the cases that the i/p-dependent potential term is included (indicated by Pwixif))) or not (indicated by J7n(x(^)))' 

kinetic function is quite similar at early times in both cases, which are relevant for the inflationary era and in both 
cases the kinetic function is negligible and consequently the same applies for the kinetic term. The same analysis 
can be done for the potential of the scalar field x(f)j both for the cases that the (/3-dependent term is kept in the 
potential or not. In Table lU we present the values of the scalar field x potential, with the (/3-dependent term included 
in the potential (denoted as V^) or not included (denoted as V^)- As we can see in Table HIl the behavior of the 


Time 

1 

0 

t ~ 10"^® 

t ~ 10“® 

t ~ 10“'‘ 


1.5 X 10"^® 

1.6 X 10“^° 

3.5 X 10"'^ 

7.73 X lO"’’ 

1/7 

3.46 X 10"^® 

1.9 X 10"^° 

8.8 X 10 “ 

9.2 X 10®® 


TABLE II: The scalar field’s x(t) scalar potential Vx as a function of cosmic time t, with the (/^-dependent term included in the 
potential (denoted as V™) or not included (denoted as VI/') 

potential for early times, in case that the (^-dependent term is included or not, is quite similar in both cases and also 
the potential is negligible at early times. 

As a final task, we shall investigate the behavior of the potential (jl]), for the chosen values of the parameters (j43|l . 
Since we assumed that the scalar field (/3 does not evolve in a slow-roll way, it’s dynamical evolution is governed by 
the following non slow-roll equation of motion, 

(p + iHif + = 0 , (50) 

dip 

with V{ip) given in Eq. (|4]) and notice that the difference between Eq. (ISOl) and Eq. (l25l) . is traced to the fact that 
the scalar field (j) is non-canonical, while (p is canonical. We solve numerically the dynamical evolution equation (1501) , 
and in Eig. [SJ we plotted the behavior of the potential as a function of time, by choosing the initial conditions for the 
scalar field ip to be (/3(0) = 10^° (sec“^), and p'{0) = 10“^ (sec“^) As it is obvious from Eig. [SI the scalar potential 
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1.4x10“- : 

1.2x10“- : 

•T 1.0x10“- : 

: 

3 8.0x10^®-_: 

6.0x10^®- - 

> 

4.0x10^®- : 

2.0x10-'^^- : 

n t, .. 

0 2.xl0-^^ 4.xl0-^® 6.X10-5® 8.X10--''® l.xlO'^^ 

t(sec) 



FIG. 5: The time dependence of the scalar field potential V{ip{t)), with ci = 10^^, C 2 = C 3 = 10^°, C 4 = 10“^'^, 

ts = 10“^®sec,a = 4/3 and for the initial conditions ip(0) ~ 10^° (sec~^), 'p'(Q) — 10~^ (sec~^), including the p-dependent 
term. 


dominates at early times and also remains quite large until t ~ 10“^sec, from which point it rapidly decays, and 
therefore the second scalar field dominates the cosmological evolution. In order to see how rapid is the decay of the 
potential V{(j){t)), in Table Hill we present some values of the scalar potential V{(j){t)), for various time points. Notice 
that at t = 10“*^'^sec, the potential decays quite rapidly and it is practically zero at t = 0. An important remark is in 


Time 

t ~ 

t ~ 10”^ 

t ~ 

t ~ lO"®--*® 

b 

H- ' 

0 

1 

0 

-4 

V{p{t)) 

7.4 X 10®® 

8.5 X 

2 X 

516.3 

4.5 X 10"®® 


TABLE 111: The scalar potential V{p{t)) as a function of cosmic time t, for the initial conditions p{0) — sec ^, y>^(0) ~ 10 ^ 
sec“^, and also with ci = 10^^, C2 = 10“^°, C3 = 10^'’, C4 = 10“'^°, ta = 10“®®sec,Q; = 4/3. 

order: The numerical analysis for the scalar field p evolution, showed that the scalar field p remains quite large for a 
wide range of time, therefore in the previous analysis of the scalar field evolution x(t), the 93 -dependent term can be 
safely neglected, since it is negligible {p —>■ 00 ). The behavior of the scalar field p{t) as a function of the cosmic time 
can be seen in Fig. [6l Finally, we have to note that the results are quite model dependent and also strongly depend 



FIG. 6 : The time dependence of the scalar field p{t), with ci = 10 ^^, C 2 = 10 C 3 = 10 ^°, C 4 = 10 tg = 10 ^®sec,a = 4/3 
and for the initial conditions 93 ( 0 ) ~ 10^° (sec p'(0) ~ 10 (sec~^). 

on the initial conditions and the values of the parameters. However, the focus in this paper was mainly to describe 
qualitatively the dynamical cosmological evolution and also to indicate that the system is dynamically unstable. This 
is the subject of the next section. 


III. INSTABILITY INDUCED BY TYPE IV SINGULARITY 

Having described the qualitative behavior of the model, we come to the crucial part of our work, the instability 
of the system. As we demonstrated in the previous section, the two scalar field model with action (fT71) is described 
by the canonical scalar field action of Eq. (|4^ . before and during the inflationary era, since the contribution of the 
second scalar field x is negligible during this era. However, the dynamical cosmological evolution of the system is 
described by the Hubble rate (1321) . In practice, the parameters are chosen in such a way so that the term responsible 
for the singularity ^ (—gives a very small contribution to the Hubble rate, which is mainly driven by the first 
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term in (1321) . Since one of our main assumptions in this paper is that both scalar fields do not respect the slow-roll 
evolution requirements, we shall study the Hubble slow-roll parameters of Eq. m, for tliG Hubblo rOito {|32J) . .A direct 
computation for the parameter ch yields, 


+^ 4 (-^ + 4 ) 
(s;^+C4(-i + ts)“) 


(51) 


while the parameter rjH is equal to, 


2 + ^4(-i + is)“) (- {c2+clt)^ ~ + is)"^+“a) 

By looking Eq. dSU, it is obvious that at early times, the contribution of the terms C 4 (—and C 4 {—t + ts)°‘ 
is negligible, compared to the other terms. We can also verify this numerically, since for t = 10“^°sec, the term 
C 4 {—t + is of the order ^ 10“^^, while the other term in the numerator of (ICTl . namely (^c 2 ‘+clt)^ 

order 10®^. Notice that we used the values for the parameters defined in Eq. (H51) . So practically we should neglect 
the aforementioned two terms at early times, and consequently the Hubble slow-roll parameter e can be approximated 
by, 


en — 


C1C3 

(C2+C3ty-‘ 

(^y 

\ C2+C3t J 


C 3 

Cl 


(53) 


and by taking into account the identifications we made in Eq. (1421) . we obtain that the early time approximation of 
the Hubble slow-roll parameter ch for the Hubble rate ([32l) is €h — |, which is identical to the one we obtained in Eq. 
(nsj. Hence, this result also verifies our claim that the contribution of the slow-roll scalar y at early times is practically 
zero. Notice also that at t = tg, the Hubble slow-roll parameter is exactly e = -, without any approximation, since 
both the terms 04 )—t -I- ts)“^+“a and C 4 (—f -I- ^ 3 )“? vanish, for the case of a Type IV singularity, which occurs for 
a > 1. Hence, although the Hubble slow-roll parameter e is small and very well defined for a Type IV singularity, the 
same issue of the graceful exit from inflation persists, that is, inflation occurs, but never ends, the slow-roll parameter 
is constant. 


1. Singular second Hubble slow-roll parameter and implications-Discussion on the results 

Let us now investigate what happens with the second Hubble slow-roll parameter tjh, for the Hubble rate (1321) . 
As we can see in Eq. (l5^ . in the case of a Type IV singularity at t = ts, and when I < a < 2, the numerator of 
rjH becomes divergent. This is very important, since this is an indication of instability in the exponential model. 
Moreover, the term responsible for the singularity, that is, C 4 {—t ts)“^+“(—1 -|- a)a, is negligible for early cosmic 

times with t > tg and t < ts, and only becomes important at exactly the singularity point t = tg. The same applies 
for the corresponding term in the denominator, so at early times and for t > tg and for t < tg, the second Hubble 
slow-roll parameter becomes, r]H — but at t = tg, rjH —^ 00 . Recall that the second Hubble slow-roll parameter 
measures how long the inflationary era is. 

So we have obtained a physically appealing resulting picture, according to which, at early times, the two field scalar 
model evolution is effectively identical to the one caused by a single scalar field with potential (|1]) , and with slow-roll 
parameters that can be compatible with observations. In addition, the inflationary era ends violently at the point 
t = tg where a Type IV singularity occurs, at which point the second slow-roll parameter rjH becomes divergent. Since 
the point tg can be arbitrarily chosen, this means that the inflationary evolution is unstable owing to this infinite 
singularity on the second Hubble slow-roll parameters at every point we may wish. 

The critical question is, what does this instability indicates? Could it indicate somehow that this could be a 
possible mechanism for a graceful exit from inflation in this exponential scalar model? The quantitative answer to 
these questions is not a trivial task however, so we discuss here qualitatively, the implications of the instability and 
the possibility to connect the instability to the graceful exit solution. 

This instability mechanism we propose for the ending of the inflationary era for the exponential model (|T]) , is quite 
different in spirit in comparison to other existing mechanisms of ending of inflation, like tachyonic instabilities [ 2 ^ 
or via trace anomaly (35j |. Note that practically the singularity we found, indicates some instability of the dynamical 
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cosmological system, but the qualitative description of this dynamical parameter rjH is that it measures when inflation 
ends, and the singularity a,t t = tg, definitely indicates strongly that the inflationary evolution is abruptly interrupted. 
So by having in mind that the only source of singularity is the second Hubble slow-roll index, it is rather tempting 
to ask if this instability could be actually indicate termination of the inflationary era. As we now demonstrate, this 
may be indeed occur. 

In order to achieve this, we shall use the so-called slow-roll expansion, which is a generalization of the slow-roll 
approximation [^. In fact, the slow-roll expansion is a more accurate approach towards finding the inflationary 
solution, in comparison to the standard slow-roll approach. The slow-roll expansion was introduced in [^ . and it is 
worth recalling in brief its basic features (for details see [l^). Consider an inflationary solution that gives rise to a 
specific Hubble rate, H{(p), where is a canonical scalar field. Note that we shall for the moment adopt the notation 
of Ref. [13, and we shall express all the quantities to be used as functions of the canonical scalar, and in the end we 
express all the parameters as functions of the cosmic time, implicitly through the variables en and rjH- In the context 
of the Hubble slow-roll expansion, an attractor inflationary solution is assumed, to which all the solutions that can 
be obtained from the potential slow-roll approximation, tend asymptotically. Note that the attractor inflationary 
solution is very important for the analysis that follows. Then, the single scalar field Friedmann equation reads. 




(54) 


where en is the Hubble slow-roll parameter appearing in (TT^ . but expressed as a function of the scalar field ip. 
By using the binomial theorem, it is possible to obtain the following perturbative expansion [^ for the Friedmann 
equation (|54|l . 


(v) (l + ev - ^4- + ^<^vVv (55) 

+ f f 4 --).- + + o.). 

where we kept terms up to fourth order, and the parameters ev,VViiv are defined in terms of the Hubble slow-roll 
parameters m, as follows. 


Tiv = {3- en) ^ (Se// -k Srjn - Vh - ^h) > (56) 

■Cy = (3 — sh) + 9Clf ~ ^^hiIh ~ ~ + ^Vh^h + Vho'h'^ ■ 

and also and an are given as functions of en and rin below. 


& = - Vh) - y■ (57) 

Note that the prime denotes differentiation with respect to the scalar field (p and also it is assumed that the solution 
H{t) is the attractor of the theory. It is conceivable that the slow-roll approximation takes into account only the 
lowest order terms, while the expansion (1551) provides a better approximation to the inflationary attractor H(t), 
which continues to describe the inflationary era, until the slow-roll expansion breaks down. As was also stressed by 
the authors of [t^ . the perturbative slow-roll expansion approximation breaks when the slow-roll parameters take 
large values, or if a singularity occurs. Using their example, for a potential of the form ~ the expansion contains 
inverse powers of p, and the inflationary era may end when the perturbative expansion is not valid anymore. This 
can occur definitely at the minimum of the potential ai p = Q, where singularities appear in the expansion, but note 
that exit from inflation occurs earlier, owing to the fact that the expansion contains inverse powers of the scalar field 

The case we studied in this paper is perfectly described by what we just presented, that is, the exponential potential 
p drives the inflationary expansion until at some point a singularity occurs, so the solution is finally interrupted and 
the new attractor must be found. Practically this means that the inflationary attractor solution may be perfectly 
described by the perturbative slow-roll expansion solution, up to arbitrary order, until the perturbative expansion is 
violated, and the solutions exit from inflation. It is conceivable that the qualitative picture of our case goes as follows: 
The dynamical inflationary evolution is governed by the solution of the exponential scalar field, since the Type IV 
singular part of the Hubble rate (contribution of the second scalar), contributes very little to the dynamics. However, 
when the Type IV singularity time ts is reached, the second Hubble slow-roll parameter blows up and the expansion 
is not valid anymore, so this could indicate that inflation might end there. 
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Recall that the slow-roll expansion is constructed on the basic assumption that an attractor solution exists, before 
the exit from inflation. So as a final task, we shall prove that the solution H{t), given in Eq. (15^ . can act as an 
attractor of the theory. Assuming that the attractor is, 

Ho{t) = —+ C4 (<s - (58) 

and that the dynamics is governed by the canonical scalar field, we linearly perturb the system H = Hq + SH, so at 
first order we get 


- {HofSH = A9 ,t:k^Hq5H , (59) 

with the solution to the differential equation being of the following form, 

5H = 6H{U)e^v{-l2K^ (60) 

where dH{ti), the value of the perturbation at some early time ti. By choosing the values of the parameter we used 
in the previous sections, and also by choosing 5H{ti) to be very small, like for example 5H{ti) = (sec“^), with 

ti = 10“^^sec, we can solve numerically the integral in Eq. (I60p . In Eig. [71 we plot the behavior of the perturbation 
5H as a function of time, where can see that the perturbation vanishes eventually, in an exponential way. Note 
however that this may depend on the initial conditions at the initial time ti, but the exponential decay is the resulting 
picture again. Hence, the solution (I55)) . is the inflationary attractor of the theory, which is identical to the analytical 
expression that can result from the slow-roll expansion. However, at the Type IV singularity, the slow-roll expansion 
breaks down, and therefore this could indicate that the solution (1551) is no longer the inflationary attractor. Recall 
again that there is a clear distinction between the slow-roll approximation and the slow-roll expansion. In the present 
case, the slow-roll expansion breaks down at second order. 



FIG. 7: The time dependence of the perturbation 8H{i), with ci = 10^^, C2 = 10 cs = 10^°, C4 = 10 ta = 10 ®®sec,a = 
4/3 and for the initial conditions 577(10“^®) = 10“^^° (sec~^) . 

Before closing this section, we need to note that if the singularity was Type HI or even Type I, the Hubble slow-roll 
index would be singular at 7 = tg, so practically we wouldn’t be able to even discuss about inflation, since the 
term ^ (7 — 7s)““^, would make the Hubble slow-roll parameter eu, quite large. In order to have a taste of how large, 
assume that a = —2 (the Type I case according to our classification) then for 7 ^ 10“^°sec, this term would be of 
the order 10^°, which is not negligible. Similar results hold true for the case of the Type HI, and for the Type H 
singularity, but a full analysis is deferred to a future work, since this would require to change all the numerical values 
we assumed for the parameters and also we would have to change the initial conditions for the scalar fields. 

In the next section with shall study the stability of the cosmological system of the two scalar helds when this is 
viewed as a dynamical system. 


A. Cosmological dynamical system Stability-Possible graceful exit from inflation via instabilities at the cost 

of phantom inflation 

In this section we shall study the stability of the cosmological evolution caused by the two scalar fields, viewed 
as a dynamical system. Particularly, we shall be interested in the stability of the solution (1701) . which results from 
the reconstruction method we used for two scalar fields. The cosmological evolution is given in terms of the Hubble 
rate dSH), and it is governed by the FRW equations of Eq. (fT51) . In order to better study the system of cosmological 
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equations as a dynamical system, we shall introduce the quantities and X, which are related to the scalar 

fields (j) and y as follows, 


X^ = 0, X^ = x, x = 
and in view of Eq. (EB, the FRW equations (fTB can be written as, 


x) 

H 


(61) 


( 


dx^ w«p)iXi-l 
dX 

dX 277 (x)X 

dX _ (1 - X2) 

dX “ X^ + X^ ^ m x^ + x^ 


2uj{4>)H 

v'ix) (x^ - 1) 


- 3 (X^ - X) 

- 3 (X^ - X) 

ff X^X^ + 1 - Y(X^ + X^) 


Notice that the solution of Eq. (EOl), can be expressed in terms of the variables dm), as follows, 

X^ = 1, X^ = 1, X = 1. 

By performing the following linear perturbations of the dynamical variables dSB, 

X^ = 1 + SX^, X^ = 1 + SX^, V = 1 + SV, 

we obtain the following dynamical system which describes the evolution of the perturbations, namely. 


( \ 

( 


0 

3 

= 


0 

'n ix) o o 

Hv(x) 

V ) 

[ 

0 

0 

-3- 


d 

dX 


The eigenvalues corresponding to the matrix appearing in the dynamical system above, are equal to. 



(62) 


(63) 


(64) 


(65) 


= — 


/(-(>) 


Hence, by substituting the Hubble rate 


-3. 


^ ^ _ v'ix) 
^ Hvix) 


- 3 , My - Jp 


( 66 ) 


we obtain. 


= -3 - 


M^ = -3- 


My — —3 “h 


2c3 


(c2+C3t 1 '=4( t + tg) ^ 

— 1 + a 

(C2 + C3(/)) 

(c2+C3t 1 C4( t + tg) ^ 
(C2+C3t)" + ^ + is) 

1 {ts - x) 




(67) 


From the theory of dynamical systems stability 37] it is known that when all the eigenvalues are negative, the 
dynamical system is stable, so the evolution of perturbations (j64p is stable. When one eigenvalue is positive and the 
other two negative, then the dynamical system is unstable, when two are positive and one negative, the system is 
stable, and when three eigenvalues are positive, then the dynamical system is unstable. Finally when one or more 
eigenvalues are zero, then we have a saddle point in the dynamical evolution of perturbations. As it is clear from Eq. 
(EH), for the choice of the parameters as in Eq. (H51) both and My are negative, but for the eigenvalue M^, it is 
always positive, if the value of x(t) is larger from tg. This however strongly depends on the initial conditions for the 
field x{t)- For the initial conditions we chose (Hll) in the previous sections, this is always true since x(0) = 10“^®, and 
the scalar field xit) increases with time, as we evinced in a previous section. Therefore, in our case the dynamical 
system always develops an instability, with this feature probably indicating that the inflationary solution driven by 
the scalar field ip at early times is not the final attractor of the system, at least for the initial conditions we chose. 
However, the full study of the cosmological dynamical system exceeds the purposes of this paper, since our aim is to 
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qualitatively stress the fact that, the Type IV singularity in some way makes the system unstable and the inflationary 
solution of the system at early times is not the final attractor of the system. 

Finally, we should note that if the initial conditions of the scalar field x(t), were chosen to be smaller than the value 
of ts, or if in our case ts is chosen to be larger than 10“^^, so that x(0) < ts, then it is possible that, as x increases, at 
some point ti corresponding to the crossing point tg = x(^i)) eigenvalue becomes infinite. This infinity may 
also indicate an instability of the dynamical system, thus indicating an abrupt change of the final attractor solution 
of the scalar field (p. From a mathematical point of view, this dynamical system analysis is a difficult issue and should 
be scrutinized in order to come to rigid conclusions. From this brief study we keep only the result that the dynamical 
system of perturbations of the cosmological parameters and V, is unstable. 


IV. LATE-TIME EVOLUTION OF THE TWO SCALAR FIELD MODEL 


In the previous sections we thoroughly studied the early time behavior of the two scalar field model, but we did not 
address the late-time evolution study of the model. Here we study the late-time behavior of the two scalar model, in 
terms of the effective equation of state (EoS) of the model. Before getting into the detailed analysis, let us recall the 
essential features of the two scalar field model we used, assuming that the parameters are as these appear in Eq. (1431) . 
Eirstly, we assumed that the Type IV singularity occurs at early times, specifically at t and secondly, the two 
scalar fields are required not to obey the slow-roll conditions. As we evinced, for the choice of the variables as in Eq. 
(H51) and also by assuming that the initial conditions of the scalar field x are taken as in Eq. (H41) . we demonstrated 
that the scalar field x has a negligible contribution to the cosmological evolution of the model at early times, but 
dominates the cosmological evolution at late times. As we now show, this behavior can also be seen by studying the 
EoS of the two scalars system. The EoS parameter Wef / for the two scalar model is given by. 


P , 2H 
WeS =-=-!- 


3Fr2 ■ 


By using the form of the Hubble rate given in Eq. (1321) . the EoS reads. 


( 68 ) 


-(C2 + C3t)^ -C4i-t + tg) 

The above relation dMl) will be the starting point of our analysis. We shall investigate the behavior of the EoS at 
early and at late times. We start off with early time, in which case, the contribution of the terms ^ C 4 {—t + ts)~^^°^a 
and Ci^—t + tg)°‘ appearing in the numerator and the denominator of the fraction in Eq. (I69L is negligible compared 
to the other terms. Also, when t ~ ts) these are identically zero, for the value of the parameter a we chose. Therefore, 
by neglecting these, the EoS at early times reads, 


(69) 


Weff = -I- 


Weff ~ -1 -I- — , (70) 

OCl 

which clearly describes quintessential acceleration, different from de Sitter though, since in our case ^ (see Eq. 
(HU)) and we took q = 100, so it cannot be considered nearly de Sitter, but purely quintessential acceleration. Hence 
the inflationary era is a quintessence era, without the possibility of evolving into a phantom evolution. With regards 
to the late-time era, it is conceivable that as t grows larger, approaching the present time era tp ~ 4.25 x lO^^sec, the 
terms ^ CA{—t + ts)~^^°‘a and C 4 (—t-l-ts)“ will start to dominate the evolution, therefore the EoS at late-time reads. 


2{-t + ta) ^ “a 

^eff — 1 5 

OC 4 

(71) 

and since a = 4/3 and also tg <C t, then. 


(-t + t,)“-i ~ (-t)“-^ = 

(72) 

Consequently, the EoS (|7T|) becomes at late times, 


OC 4 

(73) 
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The EoS (I73p describes nearly phantom acceleration, with nearly meaning that the acceleration is slightly phantom. 
Indeed, since we chose C 4 = 10“"^°, and also owing to the fact that we consider late times, when t is of the order of the 
present time, that is t ~ 4 x lO^^sec, then the EoS is approximately equal to WeS — —1.06545. Note that observations 
indicate that the EoS at present time or in near fnture may crossed the phantom divide [s^, so this result is quite 
interesting. 

In conclusion, in the context of the two scalar field model that realizes a Type IV singularity at early times, we 
were able to describe early and late-time acceleration in a unified way, by using the same theoretical framework. One 
interesting scenario that we did not address here is to study what is the impact on the two scalar field model we used, 
if the Type IV singularity occurs at late-time. However, this would require to change completely the constraints we 
used, so we hope to address this in the future. 


V. TYPE IV SINGULARITY AS A SOURCE OF INSTABILITY FOR OTHER SCALAR 
MODELS-UNIVERSAL DESCRIPTION FOR NON-SLOW-ROLL EVOLUTION 

From all the previous sections it is conceivable that the instability in scalar models caused by a Type IV singularity, 
can be achieved if the slow-roll conditions for the scalar fields are not taken into account. Then, it is possible that the 
Hubble slow-roll parameters are small during the inflationary era, but these blow-up at the Type IV singularity, where 
inflation is abruptly interrupted. So by saying that the slow-roll condition is not taken into account, it is meant that 
the indicators for the inflationary dynamics are not the potential slow-roll parameters but the Hubble slow-roll 
parameters In this section we briefly investigate if this instability mechanism we proposed in this paper, can 
also generate have implications on other scalar models if these are put into proper context. We shall study two quite 
popular scalar models, a large field inflation model [13, [IJ and a modified inflation model [H, [i3|- We start off 
with the large field inflation model, which we studied in detail in [T3| and we were able to associate it to a Type IV 
singular evolution, by using a single scalar field and not two. Let us recall the essential information for the Type 
IV singular evolution for the large held inflation and for details the reader is referred to [T3|. The Hubble rate that 
describes the Type IV singular evolution is equal to, 

H{t) = fo {ts - tr , (74) 

so for a > 1 a Type IV singularity occurs. Consider the single scalar field action, 

S = J - V{cj,) + L„,atter| ■ (75) 

where the function stands for the kinetic function and V{(p) is the scalar potential of the scalar field (j). Using 
Eq. (l75ll . the effective energy density and the effective pressure of the non-canonical scalar field (f) can be written, 

P = + ^(</>), P= -y{(/})■ (76) 

Then in view of Eq. (1751) the scalar potential U((/)) and the kinetic function w((/)) can be expressed in terms of the 
Hubble rate as follows, 

= Vicj,) = ^ (SH^ + h) . (77) 

tv tv \ / 

The single scalar field reconstruction method is based on the crucial assumption that the kinetic function and 
the scalar potential V{4>), can be expressed as functions of a single function f{(j)), as follows, 

^(<(>) = = 4 ( 3 /(</')" + /'(<^)) ■ ( 78 ) 

tv tv 

For further details on this method, see [1, [i3|- Consequently, the FRW equations for the single scalar read [sl. [l3|- 

^ = t, H = f{t). (79) 

The kinetic function and scalar potential for the scalar (p, in the case that the Hubble rate is given as in ([74]), near 
t = tg, are equal to. 


ujip) = (ts - (pT ^ 


v(<p)--^(ts-pr-^ 

tv 


(80) 



18 


By canonically transforming the single scalar field (j) (see also appendix), we can have the potential in the following 
form, 


Viv) 


afo f K(a + 1 ) 1 

1 2 V 2 a/o 


(81) 


where the field is canonical. The potential (ISTl) . is known to describe large field inflation Bm. The observational 
implications of the potential are quite appealing as we evinced in [l^, in the case that the singularity occurs 
at the end or after inflation, but in the context of a slow-roll evolution. Now consider that we abandon the slow-roll 
evolution requirement and in addition assume that the singularity can occur at an arbitrary time instance. The 
Hubble slow-roll parameters for the Hubble rate (l7^ . are equal to. 


{-t + ts) ^ “a {-t + ts) ^ “(-1-ka) 

=-A-' =- w<, -■ 

Therefore by choosing /o to be enormously large, larger than (— t -|- 4)“^““, that is, /o » (— t -I- for t ~ 4, 

then both the Hubble slow-roll parameters en and rjH can be significantly smaller than unity, eHi'HH -C 1, in the 
neighborhood of t ^ ts- But for a Type IV singularity {a > 1), both the Hubble slow-roll parameters blow up to 
infinity at t = ts, where the singularity occurs. Therefore, even in the non-slow-roll approximation, this could be an 
indication that inflation is interrupted severely due to the existence of an singular instability. Hence, our mechanism 
can work for other potentials too, apart from the exponential we explored in the previous sections. Notice that /o is 
a free parameter, so it can be chosen appropriately in order concordance with Planck data can be achieved, at least 
in the neighborhood of the Type IV singularity where inflation can be chosen to end. 

Before closing, we demonstrate in brief that the same mechanism can crea te g rateful exit from inflation for nearly 
potentials which can generate a singular evolution, which we studied in [iq . As we evinced in [l^, the Hubble 
rate for the Type singularity generating nearly potentials reads. 


H{t) = 


Cl 


+ C 4 + C 5 (tg — t) 


Cl Cst 

and consequently, the first Hubble slow-roll parameter ch is equal to, 

- (,,TcltP -C5i-t + tg)-^+^a 


Ch = 


(c4 + + C^i-t + 


2 ’ 


while the parameter r]H is equal to. 


riH = 


(,^+C5(-t + t.)-^+“(-l + a)c 


(^4 + + C5{-t + ts)“) (- (C2+C3t)^ “ + tg) l+“a) 


(83) 


(84) 


(85) 


Hence, for 1 < a < 2, the second slow-roll parameter blows up at the Type IV singularity. Therefore, in this case too, 
the Type IV singularity induced instability can work in the same way as in the exponential model we studied in the 
previous sections. Note that the model studied in and the one we studied in this paper have many qualitative 
similarities, except that in [l^ we used the slow-roll approximation for the canonical scalar p. 


VI. CONCLUDING REMARKS 

In this paper we included a Type IV singularity in the dynamical evolution of a quite well studied exponential scalar 
model, which is known to be compatible with the Planck data, when the spectral index of primordial perturbations is 
considered. In order to successfully incorporate the Type IV singularity we used two scalar fields, and by employing 
very well known scalar reconstruction techniques, we managed to produce the Type IV singular behavior appearing in 
the Hubble rate of the model. The scalar exponential model we studied however, is known to suffer from the graceful 
exit from inflation problem, but in the context of our model, we provided an indication on how inflation might end in 
this model. However, this issue has to be thoroughly addressed, a task which is beyond the scopes of this introductory 
paper on these issues. 

In order to generate an instability, we assumed that both the scalar fields, one being canonical and one non- 
canonical, do not satisfy the slow-roll requirements. Also, the values of the parameters and the initial conditions on 
the fields are chosen in such a way, so that the model exhibits the following features: 
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• The Type IV is assumed to occur at early times. In the context of our model, this will indicate the end of 
inflation in the model. 

• At early times, the canonical scalar, which is the one that corresponds to the exponential scalar model, dominates 
the cosmological evolution, with the scalar potential and kinetic term of the non-canonical field being neglected. 
Special emphasis was given on the inflationary era. 

• The non-canonical scalar never becomes a phantom scalar, so singular phantom evolution is avoided. 

As a side effect of our model, the late-time evolution is governed by the non-canonical scalar, as we analyzed in detail. 
In fact, in the context of our model, late-time acceleration and early-time acceleration can be described in a unified 
way, as we evinced. 

For our two scalar field model, since the slow-roll condition is not respected by both scalar fields, the indicators that 
determine the inflationary process, if it occurs and when it ends, are the Hubble slow-roll parameters. By calculating 
these for a Type IV singularity, we obtained the quite appealing result that at the Type IV singularity, the second 
Hubble slow-roll parameter, which is known to indicate the end of inflation, diverges at the Type IV singularity. 
As we claimed, this might be an indication that inflation ends at the point of the Type IV singularity. Away from 
the singularity however, and always referring to early times, the Hubble slow-roll indices are approximately (and in 
some cases exactly) equal to the ones corresponding to the exponential scalar model solely. This in some way also 
validates our claims about the evolution of the two scalar field system. In addition, we supported numerically and 
verified in a semi-quantitative way our theoretical claims. We need to note though that in order to obtain some 
quantitative results, the analysis must be much more rigid and also someone should examine the dynamical system 
of the cosmological equations exactly. Our aim was to indicate in a qualitative and semi-quantitative way, the fact 
that a Type IV singularity might be responsible for the graceful exit from inflation, to certain cosmological models. 
Finally, we studied analytically the stability of the cosmological equations, when these are viewed as a dynamical 
system, by examining the dynamical system of linear perturbations of certain parameters, related to some physical 
quantities. As we demonstrated, the system is unstable, with this instability possibly indicating the fact that the 
initial solution corresponding to the canonical scalar field is not an attractor of the dynamical system. Indeed, this 
was the case, since the solution lasts up to the point that the Type IV singularity occurs, with the duration of the 
solution being quantified by the second Hubble slow-roll parameter. 

We need to note that the effects of the Type IV singularity on the cosmological evolution of the model are somehow 
unique since only the second slow-roll parameter diverges and not the first. One should however examine the rest 
of finite time singularities, with the most appealing, from a phenomenological point of view, being the Type H. 
Since finite time singularities frequently occur in cosmology, the phenomenological implications of these have to be 
understood in detail, especially the effects of the mildest ones, like the Type IV and Type H. These could be the link 
between a quantum theory of cosmology and the classical cosm olog ical theory, with the most appealing candidate of 
a quantum cosmology theory being Loop Quantum Cosmology [^. Finally, it is worth to study the evolution of the 
model we studied in this paper, in the context of theories with generalized equation of state. For a recent study see 

Ill- 
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Appendix 


Here we discuss in brief why the incorporation of the Type IV singularity in the cosmological evolution of a single 
scalar field is very difficult, not to say impossible. Consider the following generic non-canonical scalar field action, 

S = I - ^(<^) I ■ (86) 

with the function being the kinetic function and the scalar potential is V{(t)). For the flat FRW background of 
Eq. ([1]), the energy density and the corresponding pressure are equal to, 

P = + ^(0), P= ^w((())(^^ -V{(()). 


(87) 
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Therefore, the potential V{4i) and the kinetic function a;((^) can be expressed in terms of the Hubble rate in the 
following way, 


^ _ 2 _h , V{(j>) = 4 (si?' + h) ■ ( 88 ) 

By making the transformation (I37|) . we can rewrite the action of Eq. (I86p . in terms of the canonical scalar field (/?, 
since the kinetic term becomes. 


- = -d^ipd^ip. (89) 

and finally the action of Eq. (l86l) reads, 

S = j - ^(<P) I • (90) 

The technique that simplifies to a g reat extent the incorporation of finite time singularities, is the scalar-reconstruction 
technique developed in Refs, [l^ [l^, which is based on the assumption that the kinetic function and the scalar 
potential can be written in terms of a function /(<()), in the following way, 

^{(t^) = — y{4') = , ( 91 ) 

Kj Kj 

Therefore, if we assume that no other matter fluids are present in the cosmological equations, the field equations (18711 
take the following form. 


^ = t, H = (92) 

As in the inflation case, which we studied in [l^ . our inability to incorporate the finite time singularity for the 
potential o, is traced on the fact that the potential has an exponential functional dependence with respect to the 
scalar field ip. In the case of a Type IV singularity, a quite general form of the function f{t) that can describe the 
Type IV singular evolution, has the form. 


f{t) = cit-tsf + fi{t) (93) 

with fi (t) being the non-canonical scalar-tensor reconstruction function that produces the scalar potential (jH . Equiv¬ 
alently, f{t) can be equal to. 


fit) = f{t) = e^* (94) 

with c being an arbitrary constant. Eor both the cases (IMl) and (IMll . it is very difficult to incorporate the Type IV 
singularity for the following reasons: 

• The function /(<) must be chosen in such a way, so that it is possible to solve the integral ip = \Jf'{(j))d(j), in 
terms oi (j) = (j){p). This is not so easy to be done for the functions f(t) of Eqs. (IMl) and (IMl) . 

• In the case that an ansatz is found and the function fi4>) is found, so that (j) = 4>ip) can be explicitly solved, 
then even in this case, the resulting scalar potential becomes severely constrained and the potential (|1]) cannot 
be easily reproduced even in parts. In order to see this, suppose that (pip) is found in explicit form, then the 
corresponding scalar potential ViPip)) would be equal to, 

= 4 ( 3 /('/'('^))^ + f'iHf))) , ( 95 ) 

and this potential must be of the following form, 

ViPip)) = + ViiPip)) (96) 

which is a rather formidable task to do, since an ansatz is needed. Work is in progress though to find such an 
ansatz. 
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